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АВ5ТКАСТ 

During recent years, logarithmically periodic antennas have been 
widely used due to their frequency response characteristics, simplicity 
of design and directivity. However, their theory of operation still is 
in a development phase, and very few models have been fully analyzed. 
The present paper is an attempt to analyze the operation of a zig-zag 
model that has the property of being symmetrical, and suitable for 
operation against ground. The radiation pattern of the antenna is 
obtained for different models of current distribution, and, finally, the 
impedance characteristics and an approximate current distribution are 
obtained, using non-uniform transmission line theory. The results 
obtained show reasonable agreement with experimental data, and confirm 


conclusions drawn from physical considerations. 
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I - LOG-PERIODIC ANTENNAS - HISTORICAL BACKGROUND 


l.l - Frequency Independent Antennas and Log-Periodic Antennas 


Frequency independent antennas were originated in 1954, when V. H. 
Rumsey [14] started the idea that a structure characterized by angles, and 
without any characteristic length dimension, should behave independently 
of the frequency of operation. As a practical example, Rumsey proposed an 
equiangular spiral structure. J. D. Dyson carried out a study of such a 
structure composed of two plane sheets of metal extending from a central 
point according to this geometry, and he proved that the model not only 
had frequency independence characteristics, but also that the interesting 
part of the antenna was confined to the area inside an "active region," 
where the dimensions become of the order of half a wavelength. That is, 
the outward area is not necessary and the structure can be truncated at a 
convenient size, and will be useful within a range of frequencies whose 
lower limit is determined by the peripheral dimensions, and whose upper 
limit is determined by the spacing at the feeding terminals, near the 
center of the antenna. 

A structure like this is bidirectional and has circular polarization. 
Also, since the active region rotates along the spiral, the polarization 
rotates too. Strictly speaking, then, not all the characteristics of the 
antenna are frequency independent, but, instead, they repeat periodically 
every time the frequency is changed so as to reproduce a reference space 
condition. The input characteristics, pattern shape and field intensity 
are constant. 

One year later, in 1955, R. H. DuHamel, who was working with Rumsey 
in broadband antenna development, proposed a structure that, instead of 


being entirely defined by angles, was interrupted by discontinuities 
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Fig. 1 - DuHamel's log-periodic teeth structure 


properly located. The first model built (Fig. 1) was formed by a 
circular sheet of metal, cut along two symmetrical circular sectors, and 
with teeth placed along the radial edges. The size of the teeth and 
their distances to the center, are related by a constant ratio p. ТЕ 15 
evident that such a structure, within its size limits, will repeat its 
radiating characteristics every time the frequency changes by p; that is, 
for all the set [p £g]. À plot of these frequencies in a logarithmic 
scale will show them equally spaced by log p, and that is the origin of 
the name log-periodic given to this kind of geometry. 

The model shown in Fig. l is a bidirectional antenna with large 
half power beam width.  Bending one of the halves of the antenna toward 
the other, so as to form an acute angle, D. E. Isbell obtained a unidirec- 
tional antenna with the surprising characteristic of radiating toward the 
vertex; that is, backward in relation to the feed. Meanwhile, proceeding 
with his work in the equiangular spiral antenna, Dyson devised another 
unidirectional antenna by wrapping the two metal sheets of the planar 
spiral around a cone. Again, he obtained a backward radiating structure, 
whose characteristics could be adjusted by means of the cone angle and 
the spiral rate of growth. 

These two pioneer models brought log-periodic antennas to the field 
of practical application; and they, and their variations, have been 
extensively used in the VHF band and above. They are particularly suited 
for military applications where more often than not, a single antenna is 
required to cover a very large range of frequencies. As it might be 
expected, these antennas have moderate gains, and for fixed circuit 
applications they do not match high gain structures. However, when a trade 


between even performance along a large band, and optimum performance only 


at one frequency or narrow band is allowed, the log-periodic models offer 


the best possibilities. 


1.2 - Unidirectional Planar Structures 

Traditionally, it is with the lower frequencies that the search for 
efficient radiating systems has been less successful. Good electrical 
performance has been obtained at the expense of size, and, although 
theoretically it would be possible to build optimum systems, the costs 
involved in the ground occupied and the hardward employed make such attempts 
prohibitive. In the HF band, the best compromise between size and opera- 
tional characteristics, during the past three decades, has been the rhombic 
antenna. Today, the rhombic still is the practical structure with highest 
gain, and when used with a range of frequencies within a ratio of not more 
than 3:1, it still might be the best answer for an HF radio link. Its 
main inconveniences are the size and the variable take-off angle, but the 
latter can be used with advantage by appropriate selection of frequency 
according to ionospheric conditions. 

The log-periodic models of the type mentioned above are not, in 
general, feasible for HF long range communications. Besides being too 
bulky, they employ sheets of metal that would make the cost too high. 

Two important developments were to come, and solve these two problems. 

The first was when DuHamel and his research team verified that only 
the peripheral strips of the metal sheets are active in the radiation 
mechanism, so that, essentially, all the structures can be reduced to 
wire frames. 

The second was when Isbell [11] built the first array of log-periodic 


half-wave dipoles (Fig. 2a). Again, here a surprising fact was verified. 


10 


Fig. 2a - Isbell's transverse dipole array 
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Fig. 2b - A tapered-ladder for operation above ground 
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The structure, as all log-periodic structures so far realized, radiates 
backward; but, to obtain proper operation, every other dipole of the 
array is fed through a reversal of the feed line connections. 

This dipole array established the applicability of log-periodics in 
the HF band. The structure can be built for vertical or horizontal 
polarization, and in the latter case, by adequate choice of the slope of 
the feed line, either fixed or variable take-off angles can be obtained. 
The first models had tubular dipoles, with radius scaled to the structure 
period. Granger Associates, who might be properly recognized as the 
pioneers in large scale building of HF log-periodic antennas, have 
replaced the tubular dipoles by wire triangular elements with appreciable 
reduction in weight, making it possible to build the rotatable long range 
log-periodic. This is a steerable antenna, mounted on a pole high enough 
so that the take-off angles are as low as 10?, and light enough to be 
rotated by a small electric motor through a 360° azimuth span. 

Incidentally, one of the problems that arose with these antennas was 
the matching of the feeding line. The antenna terminals show a balanced 
input impedance of the order of 200 ohms, and, typically, a maximum VSWR 
with respect to 50 2 below 2:1 through the band of operation. Two hundred 
ohms is not the most convenient value to carry energy from a distant 
transmitter, and the problem was that of obtaining a matching device that 
would enable the use of 50 ohms standard coaxial cable. That need speeded 
up the research in the area of baluns and matching transformers with 
outstanding results. Granger Associates have supported their achievements 
in antenna production, supplying also baluns with the required broadband 


characteristics, and multicouplers that enable the simultaneous use of one 
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single antenna with a number of independent transmitters; the total power 
being the only limitation in the system plus a small separation in the 
frequencies of operation. 

Looking for more and better, the next step in log-periodics was to 
try to obtain a structure that would use a ground image with the purpose 
of reducing its size. The reversal in dipole connections excluded that 
possibility in Isbell's array, except by the use of inverting transformers 
between radiators, which does not seem to be a very practical idea. The 
more significant successful structures are the bent log-periodic zig-zag, 
the tapered ladder, and the monopole array. The bent zig-zag, developed 
by J. Greiser and P. Mayes [8], [9], will be the main subject of this 
paper. The tapered ladder, developed by Wickersham, consists of an array 
of grounded vertical radiators (Fig. 2b) capacitively coupled to an 
inclined transmission line, the slope of which is adjusted for proper feed 
of each radiator. А variation of this arrangement, that uses triangular 
wire radiators, has been commercialized by Granger Associates. The 
monopole array, developed by D. Berry and F. R. Ore [23] uses phasing 
stubs between every two radiators, scaled with the structure period. The 
stubs are horizontal and open-circuited, so that only one wire is actually 
necessary, since its image provides the other, and the radiators are 
vertical tubes, of a quarter of a wavelength, for the respective operating 
frequency. 

In general, the half-size structures do not have as good a performance 
as the transverse dipole, and the recent advances in construction 
techniques have made the latter a stronger competitor, even for frequencies 


as low as 2 MHZ. 
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1.3 - Developments in the Iheory of Log-Periodic Antennas 


The development of log-periodic structures during their first years 
was mainly the result of good intuition and extensive laboratory testing 
of a large number of proposed models. As a rule, only after a successful 
model was obtained, was a theoretical foundation to explain its perform- 
ance developed, and only relatively recently models based on sound 
theoretical investigation are being realized. 

There have been three main approaches to the problem of analysis of 
operation of LP's. The first, which we can call the classical approach, 
uses the theories developed for the components of the antennas and combines 
them to suit the particular model under consideration. The outstanding 
example of this approach is Carrell's analysis of the transverse dipole 
array [10], who described the behavior of the structure in terms of 
currents, voltages, input impedance and far-field with very close agree- 
ment with experiments. This approach usually has the limitation of being 
valid for a particular structure only. 

Another approach, that has been pursued by DuHamel [15], [16], 
Armstrong [16], Mittra, Jones [17], [18], [19], and others is to develop 
an adequate body of information and knowledge in log-periodic networks, 
on which a general theory of the antenna problem could be based. Great 
progress has been made toward this purpose, and results obtained have 
helped in a better understanding of log-periodics, and, in some cases, 
provided new ideas for future realizations. 

The third approach, followed by Mayes, Deschamps, Patton [20], Greiser 
[8], [9], Oliner, Mittra, Jones [22] and others is an application of 
concepts developed for surface and leaky wave structures, that are 


typically uniformly periodic structures, to the active region of log-periodics 
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that can be approximated as a uniformly periodic region for limited 
purposes. Since such concepts give an excellent insight into the main 
operational characteristics of these structures, we will briefly review 
some of its basic points [4], [5], [7]. 

Consider a structure, such as a waveguide, with a periodic grating 
in one face, and suppose that the waveguide carries an electromagnetic 
wave. Due to the discontinuities K the grated M there will be 
leakage of energy through them, and, also, the fields inside the guide 
will adjust to this periodic configuration. According to Floquet's 
theorem, "for a given mode of propagation at a given steady-state 
frequency, the fields at one cross section differ from those one period 


away only by a complex constant" [7]. This can be expressed as 


where L is the period of the grating and y is the propagation factor. 
Now, if the longitudinal dimension of the guide is made to coincide 


with the z axis, and we expand the equation in a Fourier series, we get 


-j (27n/L)z Fac 


Е = Е(х,у,2)е = } E (x,y)e , (1) 


а11 п 


that is, we obtain a set of n space harmonics, each one having propagation 


constant 


If we now assume that the system is lossless, we obtain 
2тп 
ل‎ 


n Bo L (2) 


where Bo is the medium phase constant. Now, we apply this result to the 
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usual description of waves bounded by guides. In general, the fields 


are described by the wave numbers and a relation of the form 


κ «οκ, 
x y 2 
When the x = O plane is a ground plane, that is the general case 


of interest for the present purpose, the y number is zero, and, with the 


assumption of no losses, we have 


β΄ s к^ А 8^ А 82 
x Ζ 


where B = K is the free space phase constant. 

If now we recognize that the x and z phase constants are given, for 
the grated guide, by a series of the form of equation (2), we obtain the 
general form of the solution. But, before, let us consider an angle 6, 
as being the angle between the z axis and the direction of a z component 
wavefront. Then we can relate the x and z wave numbers, or the phase 


factors as 
B. = K sin 6 
and 
В, = K cos ð , 


so that the direction of leaking energy, that is the direction of the 
wavefront, is given by 


B 
cos 0 - E (3) 
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Now, according to (2) we recall that B, can take any discrete value 
В (note that n ranges from all negative values to all positive values), 
and, being so, the ratio (3) can take an equal number of values, positive 
or negative. When this ratio falls in the range -1 to +l, that corresponds 
to a solution of the cos function, and also to real solutions of space 
harsonics for which there is leaking energy. In particular, when (3) takes 
the value -1, energy is fired backward, and this corresponds to the mode 
of interest in the present analysis. 

Since the amplitude of the harmonics decente with increasing order, 


we would like to have the structure operating in a low order mode. For 


the n » -1 mode we have 
2T 
E on 
апа 
сов в = -0 - 2n/L 
-1 K 4 


In log-periodic structures, always there is some sort of guide that 
carries energy from the feed point toward the radiators, and discontinuities 
placed along the guide force energy to leak away. In the transverse dipole 
the line is the feeder, a transmission line, and the grating is represented 
by the dipoles; in the teeth structure of Fig. 1, the plane sheets form 
the guide and the teeth form the grating; the bent zig-zag is a transmission 
line deformed by the radiating elements. In each of these models, and, in 
general, in each log-periodic model, we can consider three distinct 
regions: the active region comprising the two to five radiating elements 
near resonance for the operating frequency, the feeding region, between 
this and the source where radiation is small and in general the radiators 


can be considered as a capacitive loading, and the end region, that we 


1; 


assume is not reached by any significant amount of energy, since most 

of it was radiated by the active region. Since the active region is 
small, we now assume that it is a uniformly periodic structure, and then 
apply it to the above results for leaky waves. 

Perhaps a clearer picture of such an application is by means of 
dispersion diagrams. Fig. 3 shows a dispersion diagram for equation (3), 
assuming that βρ/Κ > 1. It is a well-known fact in transmission theory 
that propagation within boundaries takes place when the phase constant 
in the medium is larger than the phase constant in the boundary material, 
i.e., when the propagation velocity in the medium is lower than the 
propagation velocity in the boundary material. For metallic boundaries 
the velocity is usually taken as that of free space, and we represent 
that condition by the phase constant K. This situation is usually 
identified by the designation of slow waves and it is said that periodic 
structures have stop bands for fast waves. These fast waves that do not 
propagate within the structure are the ones that leak away and form the 
radiated field. 

In Fig. 3 we represent the fundamental and the first and second 
negative harmonics as given by equation (2), and we show the +45° slope 
lines that separate zones of slow waves from zones of fast waves. We 
see that the fundamental is a slow wave for all frequencies, and so are 
all the positive order harmonics. The negative order harmonics are 
slow waves for low frequencies, but as the frequency increases to near 
the point B, the first negative harmonic approaches the negative slope 


line. On the line we have 


cos EI = K К 
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Fig. 3 - Dispersion diagram 


that is, we start having radiation backwards; and until the positive 
slope line is reached the radiation angle rotates from 180° to 0°. 
Meanwhile, the other negative order harmonics were having similar 
behavior in the region that we can visualize to the left of the diagram. 
Since the analysis is applied to a uniformly periodic structure, 
we only have well-defined active bands, the first of which is the 
strongest, separated by stop bands (these have the opposite meaning of 
the stop bands of transmission theory). However, in log-periodic 
antennas the structure is scaled almost continuously (typical scale 
factors range from .80 to .98), so that there is a continuous active 
band, limited only by the operating frequencies of the first and the last 
of its elements. The frequency independence characteristics of the 
antenna depend on how well its discrete scaling approaches a continuous 
scaling, and the over-all antenna performance can be predicted from this 
type of analysis. Greiser and Mayes [8], [9] have used this technique 
with remarkable success, in predicting whether or not different models of 
bent zig-zag antennas would radiate, and how the radiation characteristics 
could be improved. Mittra [17] has conducted a similar analysis for the 
transverse dipole, showing by this method how the phase reversal was 
necessary to make the dispersion diagrams approach the radiation Zone. 
Again, there is a remarkable lack of dispersion diagrams developed 
for the type of structures used in log-periodics, and, as yet, the theory 
was not really adapted for them, except under the heuristic approach 
described above. We should note here that, although this application of 
leaky waves theory in a wired log-periodic structure might, at first 
Sight, seem rather displaced, there is reason to believe that the theory 


is actually valid. R. L. Bell, C. T. Elfving and R. E. Franks [21] 
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reported measurements of fields around a log-periodic structure that 

show the existence of a forward wave progressing from the source along 
the structure, and a backward wave with a 90 degrees spacial phase 
difference. The forward wave is supported by the structure; the backward 
wave is launched by it into space. That is exactly the type of leaky 
wave behavior used in the above theory. 

Ihe mathematical tools developed to deal with the various problems 
posed by electromagnetic theory, although having in common the basic 
source given by Maxwell's equations, follow different courses and tend 
to place in distinct categories, problems that'have much in common. It 
seems to be a marked trend in recent years, to reunite electromagnetic 
theory in the sense of extending concepts and ideas developed for 
particular systems to other, apparently different, applications. 

As a final remark, we note that, since log-periodics are scaled up 
from the feeder, the active elements are the ones in the active region 
together with the smaller ones. Smaller elements in all cases tend to 
increase the propagation velocity (the feeder approaches a uniform line 
toward the feed point). If we would attempt to build a log-periodic for 
forward radiation we would scale the structure in the opposite sense, 
and would count on positive order waves for the radiation process. 
However, while in the usual scaling the positive harmonics find a high 
impedance in the first part of the feeder (they are in the stop band of 
transmission theory and mismatched to the available radiators), in the 
opposite scaling the negative harmonics would still reach the end part of 


the feeder and cause a disturbing radiation from that region. 


21 


II - THE BENT LOG-PERIODIC ZIG-ZAG ANTENNA (BLPZZ) 


The BLPZZ (Fig. A-1) was developed by Greiser and Mayes using the 
concepts of dispersion diagrams. As mentioned before, this was one of 
the first attempts to obtain a LP antenna with symmetrical characteristics 
that could be fed against ground. For that reason, it seems particularly 
suited for HF applications. 

Ihe basic geometry of the antenna is derived from a zig-zag antenna 
bent along a convenient line, so as to obtain adequate radiating elements 
and proper connecting stubs. Since the structure is progressively 
elevated from ground the connecting stubs are actually radiating, but, 
for small elevation angles that can be neglected. Their essential function 
is to provide the proper phase delay between radiators, and they can be, 
and have been, replaced by other forms of delay line, or even removed, 
without destroying the basic characteristics of operation of the antenna. 

The antenna is usually defined by the radiators spreading angle ας» 
the stubs spreading angle о > the elevation angle £, the period p, and 
the length of a wire, 2b, or the distance of one element to the geometric 
origin. In Appendix A we develop the relations among these parameters and 
other useful quantities. We note that the spreading angles and the period 
are interrelated through the slope of the wires. Although each one of 
them can be independently changed, that involves a change in the dimensions 
and inclination of the wires. Also, we note that each element, for small 
elevation angles, can be compared to a rhombic antenna made of small 
wires, and with a very large apex angle. 

This antenna has the interesting characteristic of being fed in 
series, i.e., the radiators are part of the feed line. The other planar 


LP structure so far analyzed--the transverse dipole--has the radiators 
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fed in parallel, and it might be worthwhile to compare the results 
obtained. On the other side, this one, also is a structure that seems 
adequate to analyze with conventional tools. 

In the first part of this work we obtain a set of equations for the 
radiation pattern of the antenna. Then we assign to tbe wires different 
types of current distribution, and compare the results with the 
experimental data obtained by Greiser and Mayes. 

In the last part of the work we consider the antenna and its image 
as a non-uniform transmission line and we obtain an approximation for 


the input impedance and current distribution. 


2.l - The Radiation Pattern Equations of the BLPZZ 


To obtain the radiation pattern of the BLPZZ structure above ground 
we begin by considering one element at a time. We define as an element, 
one period of the structure. Since the elements are related by a constant 
size factor, that is, the period p, the general equations obtained for 
one element are valid for the other elements as well. 

In each element, we assume that the connecting stub does not radiate, 
and we split the radiating section in two straight wires. Then we 
obtain the fields due to these two wires and their images, all referred 
to an element phase center, and add them vectorially. 

We make two more assumptions. First, we assume that there is no 
mutual coupling between elements; second, we assume that the wires carry 
travelling waves. This corresponds to considering the structure and its 
image as forming a radiating transmission line properly terminated; that 
is, without standing waves. Later we will further elaborate on this 


point. 
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After obtaining the element's equations we establish an array 
expression, accounting for spacings and feed currents. 

We outline here the development of the equations for one wire, wire 
l in Fig. 4. Details, as well as the work concerning other wires, can 
be found in Appendix A. 

In general, the distant field due to the flow of current along a 


wire is given by 


Wu шр: 1 . t 4 
E = -j7 "t p sint ORE (4) 
ol 4T 0 r 


where 0', 5, 1 and r are defined in Fig. 4, and w, k, Ho and m have the 
usual significance.  u(s) is the current distribution. In general, this 


will be given by 
LI = (5) 


where a is the attenuation factor and B is the propagation factor, along 
the wire. 

For large distances the angle 0 and the distance r are essentially 
independent of the position along the wire, so that we can consider them 
as constants in regard to s, provided that we account for the phase delay 
at the far field due to the different points along the wire. 

Again, referring to Fig. 4, if we consider the distance R from the 


center of the wire, we can express r as 


r = R = s cos 0' 


and, replacing in equation (4), and defining 
1 
lg i ШС >? 
we obtain 
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Fig. 5 - Geometrical relations in the BLPZZ 
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Now we apply this equation to wire l in Fig. 5. Неге 1 = 2b, and, then 


pong eaj [oE eee -prikte ΗΝ 
e moe o o oa ORC аы 


= =е i 
E [οι ls jk [+£ acos ө) | 





Defining 
j. —4 (VR. 
Р 2 2 kR-twt ) - е j (kR+wt) 
we obtain 
jpe {5 Ἡ ΚΕ - «55 4j] 4+) - cos 8, J^ 
E, :j&F Le e ME (7) 


[= + (E 7 60 *)] 


Since we are going to use this equation with a digital computer, this 
is a convenient form. 

Now we need to convert this equation to a system of coordinates to 
be used with all the structure. We choose spherical coordinates, as 
shown in Fig. 6, and we obtain the following transformation equations (see 


Appendix A, Parts II and III): 


соз Ө! = соз (ү-є) соз 8 - sin (y-£) sin Ө сов Ø (8) 


1 
соз -ε) sin 8 + sin -ε) cos 0 cos @ 
Ее = Еб] віп! (9) 
1 
p. = g 93 (γτε) βίη Ø (10) 
n t ° t 
9. θη sin θη 


At this point it is convenient to introduce the correction to the 
"phase centers" of the element. Here the designation "phase center" is 
used to mean the reference point of the element. The actual phase center 
was not computed, although physical considerations would lead to the 
point C in Fig. 5, that was chosen. 

Again, we assume that the angles 0 and ў, as well as the distance R 
are not altered to an appreciable degree, and we only take into account 


the phase delay in the far field. The distance c from the center of 


t? 


the wire to point C is given by (see Appendix A, Part I) 


_ b sin -€ 
“1 ~ cos ὃ о 
1 
where 
"E mt L sin £ + b cos τε 
1 b sin (γ-ε) 


The phase difference in the far field is given by (see Fig. 5) 


where 


соѕ а = ¬ [sin б, соз Ө + соз δη sin 0 cos 0] , 
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Fig. 6 - Coordinate system 





Fig. / - Correction for reference point 
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so that the factor BEES in equation (6) becomes becomes επ. 


and, finally, combining equations (7), (8), (9), (10), and this correc- 


tion, we obtain 


Ё соч) pei (E sm) pea c) 
E = ) kE L e ε EE 
CM [+005 - созе!) 
x |£os (4- €) $n 84 si (X-é) cosÓ cos $] (2) 


апа 


ik + €, COS ) «чү 455 е] [xuk (ξ Bii a) 

Кр + созе | 

pu «P. e €e sin ty-€) si 
ДАЕ 0 зон Ὁ 


EE (5) 


Now combining these equations with the similar equations obtained 


for the other three wires (see Appendix B) we have, for each element, 


Е =E, +E, +E, +E (14а) 


[т] 
I 
[т] 


(14b) 
N 
E, = ) a E (15a) 


N 
E= ) aE (15b) 
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where 8 is a factor accounting for the spacing and feed currents, that 
can be split in the product of two terms, one accounting for each. 

The feed current for the ntl element is the feed current for the 
order n element delay attenuated by the length of wires 1 and 2 and the 
connecting stub. As mentioned before, the connecting stub is assumed 
to be a uniform transmission line, so that it is only responsible for an 
additional delay. 

Referring to Fig. 5, the distance along the wires from the feed 


point of element 1 to the feed point of element 2 is 


' а (= ' 
2b, ЕЕ 2b, Hs 251 Б 251 Ξ 4b, T Ab, 4 (bı д bi) 


From Appendix À we have 


on 


' _ _1 tan as 
51 


Б {5 tan aE 


so that the distance becomes 


tan as | 


Ур tan ak 


The distance to the feed point of element 3 is 


ιο as a 


Vp tan aE 


tan as 


Vp tan ak 
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and the distance to the order n element is 


tan as 


Vp tan ak 





d = 4b 
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where SUM » 1 nm Es ОЕ 1 : 
р 2 п-2 
р р 
If we subtrace 5 from S, where, 
wot ыша чанында ЧА Ыы, 
> p^ p 
we get 
-1 
M cul дур ч 
sQ p ) = 1 P те 1/р 
р 
and the distance becomes 
п-1 
d = 4b, 1 + Fan as se 
Yp tan aE р 


If we call 8 the propagation constant along the radiating wires, and 
assume that along the stubs the constant is K, as in free space, we have 


the phase angle of the feeding current given by 


2d 
8 , tan as 1 291 /p- 
4b. K К ш = IBN (16) 


Vp tan aE 


We could obtain an analogous expression for the attenuation, assuming 
that the attenuation factor was constant along the radiating structure; 
however, physical considerations imply different attenuations along the 
line, so that we obtain, for the order n element feed current an attenua- 


tion given by 


) 4 αι (17) 


where a, is the average attenuation for each element. 


Then the correction factor for current in the array expression is 
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т-і lan %s лу! 
- дың - j4 iK (E + Z ο γε бє) 
е Cz 18 


To introduce now a correction for spacing, we are going to use as 
the reference point the "phase center" of the array. Again, we do not 
know exactly where the phase center is, but we expect in that respect a 
behavior similar to other log periodic structures [10], [21]. Being so, 
the phase center should lie by the element that is nearest to resonance, 


and that we call the M element. Then, according to Fig. 8, we obtain 


= -- = " 
R. R (Ly COS ε L cose) cos © 


where cos 0" = sin 0 cos @ and 
R =R- (L сү; ) sin 0 cos Ü cos є (19) 
n M n 
In equation (6), aes becomes CREE ας αρ απ зїп о созро 


and we insert in the array expression the factor 
ο Κ( L )sin 0 cos Ø cos e (20) 


Then the total correction a, 15 


m-l ; ten xs -Yp іт со Os 
-5 4b x, -jábik (È + FE ) = +jk (La-L a) θεοοφς a 


a= et 


Finally, the array equations become 
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tan сш s in Û cos co 
N EX дье της ejf ec Coin ia 


"zl 
(22а) 
Esti) Ee sae et amet 
da š x Eg, 


(225) 


And from here we obtain the total field E by taking the real parts of E 


and δρ and combining them vectorially, so that 


-I Real (Eo) 
E (0,6) = [Real (Ег) + Ге EY — / 2” ot Ee) 
. lel LY 


where Y is an angle counted positive in the clockwise sense starting in 


the direction of the ag vector by an observer at the array phase center. 


2.2 - Current Distribution - Radiation Diagrams 


Before programming the equations obtained in section 2.1 in a 
computer, we had to decide on the current distribution that would be 
most closely approximate to the real situation. With that purpose in 


mind, we made the following considerations: 
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Note: n denotes resonant element 


Fig. 9 - Models of current distribution 
(1) As obtained in other LP networks 
(2) Assuming travelling wave attenuated at 
active region 
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(1) The antenna and its image form a radiating transmission 
line, and since experiments [8], [9] show that its input impedance 
is mainly resistive and almost constant along the operating band, 
having a small capacitive reactance probably due to input termination 
effect, it is reasonable to assume that the line is well terminated. 

(2) The elements forming the feed region must see the same 
impedance toward the end, and that is the input impedance of the 
antenna. This is so because that impedance does not essentially 
change with frequency, and a change in frequency is equivalent to 
a change of antenna. That is, when we decrease the frequency by one 
period, we move the active region one period toward the end, and 
that is equivalent to using one antenna "one period larger;" the 
impedance that was seen by the input terminals must now be seen at 
the end of the first element. But since the input impedance remains 
the same, that means the first element is a good match, and so must 
be all the elements forming the feed region. Then the current along 
the feed region must be essentially constant, except for a relatively 
small amount of attenuation. 

(3) The energy reaching the end region must be very small. The 
structure is abruptly terminated without any dissipation load, and 
it has been shown experimentally [8], [9] that the effect on the 
radiation pattern of ending the structure by a short or an open 
circuit is very small, for frequencies well above the lowest frequency 
of operation. Then, it seems reasonable to assume that the current 
dies off after the active zone, and is essentially zero at the last 


elements. 
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(4) The current distribution at the active zone is more 
problematic. The non-uniform line forming the antenna is continuous 
in the sense that its impedance [2] changes without interruption, 
and has the property of a constant product of reactance times 
susceptance. We can call this a continuously tapered line, with 
linear taper, of the type often used to match impedances. Under 
this point of view it is continuously matching itself. 

On the other side, if we think of the line as composed of series 
sections connected by stubs, then we have to consider the fact that 
some of these sections (the ones forming the active region) are near 
their natural resonance frequency and are likely to form standing waves. 

The first point of view would lead us to a heavily attenuated 
travelling wave. The second leads us to the type of current distribu- 
tion found by Carrel [10] in the dipole array, and by Mittra and 
Jones [18] for a theoretical line. A point to note here is that both 
these cases had parallel loads, while the present case is of series 
loading. Fig. 9 shows the two types of current distribution corre- 


sponding to the two points of view. 


We started by trying these two models. The results obtained can be 
summarized with reference to Fig. 10, 11, 12, that show diagrams obtained 
with the data given for the antenna referred to by BLPZZ 19c in Greiser's 
and Maye's work. This antenna operates from 500 to 2500 MHZ. Fig. 10, 
for 2000 Mcs is what can be called the typical pattern. Most computer 
runs show that kind of diagram. Half power beam widths of the order of 90° 
in azimuth and 35° in elevation. The directivity, computed by the 


approximate formula given by Kraus [6], is around 7.5 db. Also, in almost 
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Fig. 12 - Radiation pattern 
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every case the difference between the results for the two types of 
current distribution is small, although the distribution of curve 1 in 
Fig. 9 gives slightly better results, particularly in what concerns the 
front-to-back ratio. 

After an extensive set of runs was obtained, it was found that for 
certain frequencies the diagrams depart appreciably from what was 
expected. Two extreme examples are shown in Fig. 11 апа 12. The first 
is "too good," and the second is "too bad." The most likely explanation 
that we can find for that is in the way the current distribution is 
processed in the program. The program has a sample of the curves in 
Fig. 9 and assigns a value to each element. Now, that sample is always 
the same, but the real situation in the antenna varies with the frequency. 
When it happens that the wires forming one element have length of one 
quarter of a wavelength the sample fits exactly, but when that dimension 
"falls" somewhere between two elements, it is clear that we are doing a 
rough approximation, because in that case the program takes the element 
just below it, and treats that element as if it were a quarter of a wave- 
length. 

To explore this point further, we made a set of runs at closely 
spaced frequencies, covering one period of frequency values. Table I 
shows the results. For each frequency, column 2 shows the corresponding 
quarter of a wavelength, columns 3 and 4 show the sizes of the wires 
nearest to that value, columns 5 and / show the obtained directivity and 
orientation of the E vector all around the ground plane for the current 
distribution of curve l in Fig. 9, and columns 6 and 8 show the same 


results for the current given by curve 2. 
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Wire's length 





820 


855 


875 


895 


935 


970 


‚091 ‚078 
. 088 .078 
. 086 ‚078 
‚084 «078 
«080 «078 
„077 ‚067 


2092 


1992 


.078 


Current Dist | Polarization 
(1) (2) | 
DIR/ISO DIR/ISO; (1) (2) 
ο. 5.00 | -90? -90? 
1.45 7.05 | του +90° 
6.36 6.06 +90° +90° 
5724 3.9 +90° +90° 
+90°within |-90? within 
8.89 10.9 +60° +60° 
+90°remain. | +90°remain. 
directions |directions 
-90°within |-90°within 
+90° remain. |+90° remain. 
directions |directions 
Table I 


Variations in Results Along One Period of Frequency Values 
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Here it becomes clearer that the current distribution 1 results in 
a more balanced set of results than distribution 2. It is interesting 
to see the way the direction of the field vector changes along the period. 
Actually, this change of frequency is equivalent to moving the structure 
one period in space, at fixed frequency, and, in any case, we would 
expect the change in phase that has been observed in all LP's. It has 
been called the phase rotation principle--the phase of the resultant field 


changes by 360° with a change of one period, in space or frequency. 
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2.3 The BLPZZ Antenna as a Non-Uniform Transmission Line 

With the purpose of obtaining a better description of the current 
distribution in the BLPZZ, we consider the structure and its image as 
forming a non~uniform transmission line. For easier handling of the 
equations involved, we split the structure in sections, each composed of 
a radiator element and the following connecting stub. Then we apply to 
the radiator Schelkunoff's equations for slightly non-uniform lines, 
and assume that the stubs are uniform lines with characteristic impedance 
equal to their average impedance. To account for radiation, we introduce 
in the equations an attenuation factor and a phase factor obtained from 
the radiation resistance of each element and its average impedance. The 
value of radiation resistance used is one obtained by Chaney [24] and Klamm 
[25], for rhombic antennas in free space. 

According to Schelkunoff [2] a non-uniform transmission line can be 


described by equations of the form 


V(x) 


VQ (x) + νι) + V, (x) + ... (25a) 


I(x) 


Ig G0 + L, (x) Р L, (x) "ws e (25b) 


where Vo (x) and Ty (x) are the expression for a uniform line with the 
average characteristics of the line under consideration. When a line or 
a section of line is slightly non-uniform the terms of order 2 and above 
can be neglected, without large error. In particular, when the line has 
the property of having a constant ZY product, where Z is the line 
reactance and Y is the line susceptance, and Z and Y vary slowly about 


their average values, the first order term is given by 


45 


V, tx) = Vy L Box) cosh [x - A(x) sinh rx) К, Il. [Aw cosh [x -B &) sinh Γκ] (26α) 


| + М (Вози Авен + ГАН знг, -Beocosh|y] (26b) 


о 


where 


к= Ул is the propagation factor 


Κοςν Ζρ/1ρ is the average impedance of the line 


x ^ 
A(x) == Í Z (s) cosh 2Ts ds 
K 
0 0 
ια; 
B(x) == | Z (s) sinh 2Ts ds 
K 
0 0 
Z = Z — Z 
av 
Zo = 27 is the average reactance of the line 
Yo - T is the average susceptance of the line, 


and we assume that the wires are perfect conductors, perfectly isolated, 
so that there is no resistance or conductance to account for. 

The required condition of slow variation of line parameters can be 
met by considering conveniently small sections of line at a time. 


Using the hyperbolic form of the equations for a uniform line, we have 


Vo (x) Е Vo cosh [x - Kolo sinh [x (27a) 
M 
І (х) = - — sinh [x + I, cosh Ix (27b) 
0 Ky 0 
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and, then, 


Ved =V, [U+ Bi) cosh x - Aeg sih ry] - Ks Ly [A o0 coshf x &(1- Boo) si^ fx ] 
(282) 


Е -M [01+ Во) А Гх-А(х)со5Ьх]+ 1, | At) sinhl x + (3:19) cosh Px ] 
Ko (28b) 


For the line parameters we use the equations for the inductance 
and capacitance per unit length of a two wire line, in the logarithmic 
form, implying the assumption that the radius of the wires is very small 


compared with their spacing, and, so, 


Z = j bain ( =) = jel20k Ine@e ) (29) 
Ү= ј мє т/а (2) = j k/120 In ( 2) (30) 
where 
w = 2m times the frequency 
w= 47 x Πο’ Henry/m is permeability of free space 
ε Ξ 10 7369 Farad/m is permittivity of free space 


s is the spacing between wires, and 
a is the wire radius. 

The propagation factor, I, is a constant if we do not account for 
losses. If we consider the radiation resistance as a distributed loss 
resistance, with value Ra per unit length, we can use the following 
approximations described in the literature for lossy transmission lines, 


and define an attenuation per unit length. 


„ака 
аа 2K (31) 


and an average phase factor 
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where 


lo 


k is the free space phase constant, and 
M TOES 

As mentioned before, we use for the computation of radiation resist- 
ance the equations given by Chaney and Klamm for a rhombic antenna in 
free space. Those equations were developed for an antenna well terminated, 
and carrying a travelling wave. Actually, the present structure does not 
exactly have the rhombic shape, but the elevation angles considered are 
small enough to allow the assumption that they can be neglected. More 
daring is the assumption of a travelling wave. However, from the work 
reported by DuHamel on nonuniform transmission lines [16], we conclude 
that the value of Q should not be expected to be critical, and that is 
essentially what would change directly with the value of radiation resist- 
ance. 

So far we have not accounted for the mutual impedance among radiators, 
and, naturally, that is due to the large amount of computations needed, 
that would exceed the time allowed for the present work. After presenting 
the results obtained with the present assumptions we will try to do a 
rough estimate of the mutuals effects, based on the values found for 
current distribution and impedances. 

Using equations (28) for the radiators and equations (27) for the 
stubs, we start from the end of the line, where there is a short-circuit, 


and proceed toward the feed point computing the input impedances, and the 
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coefficients of the equations. Аз а block the equations can be expressed 


by 
VG) =A Vo - Ko B Lo (33a) 
o s V NEN, (33b) 
κ. ο 0 


and, for each section, we obtain A, B, C and D, and the ratio νο/1ρ: 

When we reach the feed end of the line ratio Vo/ To is its input impedance. 
Then we assign to the current the value 1.0 + 10.0 and we proceed forward, 
using the previously computed values of A, B, C and D, obtaining the 
values of current and voltage along the line, for as many points as the 
number of sections considered. From these, and using the equations, the 
complete description of the line can be obtained. 

We program these equations into the IBM 360 system, according to the 
program in Appendix C. Appendix D shows the details of the equations 
involved. 

Sample program SMIGUEL shows the results for the same antenna that 
was used in sample program AZORES (radiation pattern), for three frequencies 
within the operation band--1100, 1700 and 2300 MHZ. For each frequency 
the first run does not account for mutual impedance effects. The second 
run includes the mutual resistance that would be found between parallel 
wires with sinusoidal current distribution. 

The experimental results obtained by Greiser and Mayes [8], [9], with 
the BLPZZ fed against ground show that the input impedance has a resistive 
component varying between 90 and 250 ohms and a reactive component 
varying between +10 and -100 ohms. Since we are considering a double 


structure, we would expect the resistance values to be about doubled. 
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Although the behavior of the reactive component is more difficult to 
predict, we would expect the same orders of magnitude. 

These assertions are confirmed in the results obtained for 2300 MHz. 
Referring to SMIGUEL the column of ZINR values shows the input impedance 
seen at the feed end of each radiator, and, not only the input impedance 
is as expected (291.3-j314.8 ohms), but also, as assumed before, each 
element essentially provides a match for the next element. The column of 
CUR values shows the magnitude of the current at the feed end, at the 
vertex, and at the stub end of each radiator. We note that, for this 
frequency, element 2 is smaller than a quarter of a wavelength, and that 
element 3 is larger than that. Up to the third element the current 
magnitude values are typical of an attenuated standing wave situation. 
After element 3 the distribution is of a heavily attenuated travelling 
wave. 

As mentioned beofre, the more delicate assumption in this approach 
was to apply to this structure the radiation resistance given by Chaney 
and Klamm [24], [25] for the rhombic antenna. It appears from these 
first results that such assumption is reasonable for the large elements, 
but it is not adequate for the small elements. Also, for these elements, 
the effects of mutual impedance become increasingly important, since 
they are "electrically closer" to each other. These two sources of error 
become more important as we lower the frequency. At 1700 MHz (the 
resonance dimension falling between elements 7 and 8) the results 
obtained for the first 4 elements are inconsistent. 

A closer inspection of the results shows that the departures in 
impedance values are related to extremely high values of Q (above 20), and 


that fact even more strengthens the idea that the loading due to mutual 
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effects plays an important part in the behavior of the "feed region" of 
the antenna. 
At this point we found it useful to make some numerical analysis of 
radiation and mutual resistance values, with the following conclusions: 
(1) The radiation resistance of the "rhombic" elements increases 
steadily with size, from a few ohms at the "feed region" to 
thousands of ohms at the "end region." 
(2) The mutual resistance between parallel wires in adjacent 
elements with sinusoidal distribution of current, in phase, decreases 
steadily from hundreds of ohms at the "feed region" to essentially 


zero at the “end region." 


This enables us to obtain a better approximation, by also including 
the mutual effects, taking two parallel consecutive wires of the structure 
at a time. Again here, we neglect a few points; namely, that the parallel 
wires considered do not have the same size, the input currents are not in 
phase, the current distribution is not exactly sinusoidal, and that there 
are other interactions to be accounted for, among the non-parallel wires. 
However, before taking this step we made extensive numerical calculations 
with different arrangements of wires with different conditions of phase 
differences and current distributions, and the results obtained were all 
within the same orders of magnitude. 

Essentially, what we are doing here is to introduce a rough correction 
to lower the Q of the "feed region," and the same result could be obtained 
if we assigned to it specific (and reasonable) values. As shown in the 
referenced paper by DuHamel [16] the exact value of Q is relatively 


unimportant. 
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We note that, by this procedure 

(1) at the "feed region" we are essentially neglecting the 
travelling wave approach since the rhombic radiation resistance is 
much smaller than the mutual resistance, and 


(2) inversely, at the "end region," we are neglecting the sinus- 


oidal distribution approach. 


The results obtained with this procedure apparently prove its 
soundness. For the three cases previously considered, the input 
impedance values now fall in the expected range (439.1-j148.5, 562.1-j109.8, 
and 516.4+j19.9 ohms). The current distribution shows the same features as 
before but with stronger attentuation, making it resemble more closely an 
attenuated travelling wave configuration. 

The next step would be to obtain a more exact description of the line 
by considering a larger number of smaller sections, obtaining better 
approximations for the non-uniformity corrections, and a more detailed 
picture of the current values. That would enable the application of these 
results to the radiation equations. Unfortunately the attempts to program 
this detailed analysis failed due to unknown reasons, related to the 
difficulties encountered during the "adaptation" phase of the new IBM 360 
System, recently installed at the Naval Postgraduate School. Although we 
have received extensive help from the Computer Facility staff, and 
undoubtedly the program would eventaully work, there was not time available 
to pursue the project. 

Before closing this subject we note that, due to the method used in 
deriving the coefficients in equations (33), there is a cumulative source 


of error in the approximations used. This is so because every new value 
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is based on another value that had been obtained by an approximation, 

and only the first value is an exact value. Since we start from the end 
of the structure, the "end region" coefficients are little affected, but 
the ones at the "feed region" are likely to be farther away from exactness. 
This might be another reason for the worse results obtained with the 
smaller elements. The way to correct, or at least to minimize this error, 
would be to split the line in more sections, as was attempted. 

Finally, we note that the current distribution does not show the 
marked peak near the resonance region obtained with parallel feed type of 
structures, or circuits. On one side this seems a reasonable result due 
to the different physical lay-out of the antenna; but, on the other side, 
it remains to be explained why the radiation patterns show acceptable 
results with different current distributions. An inspection of the phase 
distribution of the current (column CANG) shows that in all cases the 
current at the feed point of any two consecutive elements is very nearly 
in phase opposition, except for the element with dimensions closest to the 
resonance dimension, and the one before that. For these, the phase 
difference is near 90°, as also happens in other LP structures. This means 
then that only the “active region" is important for radiation characteristics, 
and the two models of current distribution previously used are essentially 


equivalent to each other and to the computed values. 
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APPENDIX A 


I. Geometric Relations in the BLPZZ 
With the aid of Fig. A-1l the following relations are evident: 


L b 
on Pn 
L RT and b REI (A-1) 
P P 








— 


The relation between the length of the radiating wires and the stubs 


can be obtained in the following way: 


2b cosy 2b" cosy 
Z= Se 1 LL а Δ 
τη tan oF L, tan as (A 2) 
2b, COSY 2b _1 cosy 
x «0 —— = = ° 5 ' à 
"Να === m Lan OF 2b__,siny + 2b_siny + 4b у siny 
cosy b 
п-1 r T ] T ‚ 
tan aE ( p 1) = δ} siny (1 + р ) + 2b _1 siny 
" 
n - = 1 
tan aE E b _1 tanyY ΓΕΝ; 2b, , P tany 
Pn 2p tany tan aE 
. 4P tany tan ob T 


b — (1-p) - (1*p) tany tan aE 


Other useful relations are 


\ 


| - EXEC 
2L _1 L -1 at: L _1 tan aE tany + L. L. tan a any 


а 


= L eco fL 


Т2 (1 - >) tan aE tany 


п-1 


= ET 





[(1+p) + (p-1) tan aE tany] 
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a) Basis structure 
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b) Bent structure 


Fig A-l - Geometric Relations in the BLPZZ 


from which 


t p. boc 


δ} mmm [(1+р) - (1-р) ёап oE tany]l tan os 


and combining with (A-3) 


bai 2р lany lanxe 


Laa ер) ар) mae any] lan x5;  U-P)-0FP) lany An xg 
4 bcosy 


2b cosy 


ς I ni 
and, since L -1 TX 


(!-p) - C *p) lany lange = (I+p) lony lanes - (1-p) lam exe lanxs tan*y 


(LEÊ) (tno stones) = | (Ê (lanes tomas) lanag tanas 


lany = (4-4) 


2 lancye lanes 


Since the structure is usually defined by the parameters aE, as, ›, 
b for the last element and the number of elements, equation (A-4) gives 
the next parameter, y, and equations (A-1) and (A-2) define the rest of 
the structure. The length of the stubs can be obtained with (A-3), or 


approximately by 


с 


tan as 
b' .— ae (A-5) 
n tan aE 

Ур 

This equation is exact when as = aE, but it can be used as a good 
approximation for the typical values used. The equation is obtained by 


"optical illusion." 


inspection of Fig. A-l, with the help of a sort of 
In fact, when as = aE the structure repeats itself at Ур increments, with 
an upside down rotation. When as # aE, however, the rate of change of b' 
with as is not linear and the equation is not true. 


The relations concerning the "phase center" of the element are 


obtained from Fig. 3, by inspection 
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1 ó = I Siné + 6 COS (x-€) (A-6) 
ο. OS 


lan ο ln sn € + 6н соѕ (0+) 


A-7 
b,, sin (yté) aie 

C at 5, $m (Y-£) 
| Cos Ó, (A-8) 
Cas Š 90 e (A-9) 


II. Transformations of Coordinates 
Expressions for the angle between a wire and a direction -6'. 
Considering Fig. A-2, that is a planification of the Y=0 and Z=0 planes, 


we obtain the following relations: 


cos a sin b = cos8 


cos a cos b sin0 cos 9, Q0, = 180-0 


il 


R cosé' R cos a cos (8-b) 


li 


R [cos B cos a cos b + sin 8 cos a sin b] 


R [cos B sin 9 cos 9, + sin 8 cos 6] 


R [sin B cos 0 - cos 8 sin Ө соз 6] 
The wire in the figure corresponds to wire 1 in a BLPZZ element, and 


the angle 8 is given by (see Fig. A-1) 


B = e + 90° - ү = 90? - (ү-є) 


99 
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Fig. À-2 - Transformation of coordinates 
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Ζ΄ 


U ο 
`Р 


so that, for wire 1] 


cos O, = cos (Y¥-£) cosÓÜ -sin (y-£) Sin 0 cos 2 (А-10) 
For wire 2 we have instead of an angle 8, and angle -8, given by 
-Bos — (90 = (үе 


and using the equation for cos θη we get 


£05 0; —- Cos (f*£) cosB - sin (Y+E) Sm 0 cos Ø (A-11) 


For wire 3, we obtain (see Fig. A-3) 


cos a sin b = cos Ө 


cos a cos b sin 0 cos Q 
9 


cos? = cos a cos (8-b) 


cos a cos b cos 8 + cos a sin b sin B 


= cos 0 sin 8 + sin 0 cos 8 cos M 


and, with ó, = апа В = 90 - (ү-є) 


cos0; = cos 0 cos (ү-є) + sin 0 sin (γ-ε) cos @ (A-12) 


For wire 4, we note that the only difference is again the direction 


of the angle 8 = 90 - (ү+є), so that 


cos0, = -cos 0 cos (yte) + sin 6 sin (yte) cos Ø (A-13) 
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Fig. A-3 - Transformation of Coordinates - Wire 3 
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III. The Components of the Е Vector in Spherical Coordinates 


According to Fig. А-4, we obtain the following relations, for wire 1 


Eg = fg 60А (А-14а) 


sinA, (0 


Using the expression from spherical trigonometry 


cos a - cos b cos c 
1 sin b sin e 


we obtain 


cos A Е cos (4o- 8) - cos 0, cos 0 й sin 8 - со50 (5/0 Вс050-с058 5іпдсоѕ#) 


Sin 6, 5m 8 |-соѕ9' 5/ Ө 


sin B - sin со5°@ + cosB SING С05@ cos ff sin 30- ©05°8) +6058 5п0с050с050 


у [-со50! 570 Ji-cos*6 sin @ 


sin B sinb +t cosp cos 0 cos 


V 7- costo 


and, since 8 = 90 - (ү-є), 


соз А, 495 (0-Е) 56+ 51 (0-6) созд соз 0 


A-15 
ν 5 cos? o; ( E 


es 
We have the denomination in the form V I- cos Q, because in the field 
expression we have a sin 0' in numerator that cancels with this. Now we 


obtain, from here, sin A 
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Fig. A-4 - Components of the field vector 
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e ё enu 256 - σα 
sin'A, N I- cos* A, S /- (sin^B cos* 8» cos"B sim*8 cos - 25m BcosB 5in8cosó cosgi) _ 
Sin? o, 


(558 sm *2 9 Cos^ B cos*8 cos! gl * 2 $inBBcos B Sm 0 Со @ cosg ) 
S» $07 


that reduces to 
SwA, . Sint p- cos? ø (1-sin B) _ (- sin*B)(I-ces) _ cosg Sin Ø 
51° д) ` sin Oy sin? 6} 
and 
Sn A, - Cose sin Ø š Sin (y-£) Sin Ø (A-15b) 
om 2 56: 
The procedure is analogous for the other wires, and the results are 
for wire 2 
Е Ξ-Ε cos A (A-16a) 


0 ө' 2 


9 Е ө sin А, (A-16b) 


[т] 
li 


соѕА . Cos (£F *E£) sin @- Sin(f*&) со50 соз 
sin 6j (A-17a) 


Simia o Simi hip Sum 2 (A-17b) 
sm p 
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for wire 3 


and for wire 4 


IV. The Angle a for the Correction of Reference Point 


sin A3 


cos A 


sin A 


Ë 


sm 03 


--E 


= Eo: sin A 


n 4 


cos À 


5% (0-Е) sing 


sin A, = sin(yte) 5m £ 
sin 04 


cos Ay -~ 605 (¥-€) SinG- Sin (¥-) cos 0 cos # 


Sin δὲ 


cos A, <= ©°5(д+є) 508 * sin (y*£) cos0 cos d 


' / 
$m €, 


(A-18a) 


(A-18b) 


(A-19a) 


(A-19b) 


(A-20a) 


(A-20b) 


(A-21a) 


(A-21b) 


This is the angle between the direction of c (see Fig. 3) and the 


direction to a point in the far field. 


Since c is in the same vertical 


plane as the wires, we can use the relations developed for the conversion 


of the angle 6' to spherical coordinates, with ô in place of Вв, ара ме 


obtain the following relations: 


COS 


COS 


COS 


COS 


| 


i 


-[sin б, COS 


sin ó, cos 
1 


sin $1 COS 


sin ὃ, COS 
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COS 


COS 


COS 


COS 


sin 


Sin 


sin 


sin 


cos 0] 


cos 


COS 


COS 


0 
0 
0 


(A-22) 
(A-23) 
(A-24) 


(A-25) 





с) ме 5 d) wie 4 


Fig. A-5 - Reference point and direction angles 
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The angles б, апа ὁ, are shown in Fig. 3, and given by equations 


(A-6) and (A-7). Fig. A-5 shows the way how the o angles are defined. 


68 


APPENDIX B 


Field Equations for Wires 2, 3 and 4 
I = Wire 2 

Since wire 2 (see Fig. 5) is in series with wire 1, we need to add in 
the exponent a factor corresponding to the length of wire l - 2b - and we 
get from (7) 


ja8 — [00-0 6:)]Ь -[« «9k E - es 8; )]b 
s Le سے‎ F (B-1) 


ш куа [<+ (ЧЁ. - соғ 6:)] 


Now, using equations (A-11), (A-16), (A-17) and (A-23), we obtain 


AF ME b С. созе, | loca je - -<050,))b -| [jy - cos Ó 3 
lE 
(«+k (f -<os 4)] (B-2) 
X [cos(yré) sin 8 - $i (y*£)cos8 cosf | 


kF jk] зЁ нса созо [jk (Ê - <s 6:) - jk (# -ces 6;)]o 
Е Б ел as si 


asii (οτε) sin ZE 


-jk c š Р 
J 2 is the correction for the reference point, 


cos “a 


where the factor e 2 


being Co and cos Oo given by equations (A-9) and (A-23). 


E — Wire 3 

Here we consider the current flowing upward, but to use the same 
exponential representation for the current, we consider the distance as 
negative away from the source. This is the usual representation for trans- 
mission lines, inherent in the present assumption of a travelling wave in 


the antenna. Then 
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апа 


-jpe [5 JH $605 63) \6 [e CE + cos 0, )] b 
αι κ. ο T o шры 


Е ЕКГ зё е е a 
03 Fk J, purs fact )k(-f +605 85) (B-4) 


Now, using equations (A-12), (A-18), (A-19) and (A-24), we obtain 


ΚΙ ὁ + 63 cose] [« +) KE + c05 &;)]b -[ec (Ê +05 6)]L 


Eo, = Ek le е | 
X [cos (y-£) 5n 8 - 5in (y-£) cos 8 cos Z| (B-5) 

E [£ ° +3053 [ |5 kÈ + 05 2) -| +} К (в + CoS 8;)] b 

Eg, =~ Fk Le е  -¢ . TT х 


x Sin (J-E) sin Ø (ваш 


where е» апа a, are given by equations (A-8) and (A-24). Note that 
Cy = Cys and Cy = Cy (Tito 5). 
III. Wire 4 


This wire is in series with wire 3, so that the procedure is analogous 


to wire 2 that was in series with wire l 


-} p36 [x+jk(Z +605 e )]5 + +605 64 )Jb 
K -E E cO. νιν 


E,- Ek I smo e (B-7) 
D, ρα 4 [<+ j k (Ê +٥504 ( [ 


Now, using equations (A-13), (A-20), (A-21) and (A-25), we obtain 
TEL $64 cose] [x «9k 1 1958/)] e jk (Ë +cos 94 )]b 


Eg, =- Ek hg [x +K νεος ο) 2 


x [ cos(yté) 5m8 * $m (q*4&) cos8 cos gd | 
(5-6) 
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APPENDIX C 


Equations for Non-Uniform Line 
I. Line Parameters 
Neglecting wire resistance and conductance between wires, the line 


parameters are 


N 
li 


j 120 k 1n (s/a) 
and 


Ү 


j k/ (120 1n (s/a)) 
as given by equations. 


From here we obtain the line impedance 


2_ 


K = yZ/Y = 120 1n (s/a) = 3k 


and the propagation factor not accounting for losses 


Г = /ZY = jk 


In terms of the structure defining dimensions and angles (see Fig. 
A-la), we obtain the following equations 


For the radiators 


Z (x) j 120k (B + 1n (A + x)) for ο ος ob (0-18) 


Z(x) 


j 120k (C + ln (D - x)) for 2b « x « 4b (C-1b) 


where 


1 - 2b siny , 
А = sine 
cos (Υ-ε) 


2 cos -ε 
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2 cos +є 
n 


d 2b sin y | 20 сов (ү-є) + 2b 
= соз (ү+є) “Gae cos (ү+є) 


For the stubs 
z (x) = j 120K ( E + Iin (Hx) ν) (ς-2) 


where 


E - 1n 25іп ysin e€ 


a 
4 L L 2b tan w. = 
Jp sin Y Jp tan ας 


and the quantities involved are as defined in Fig. A-la 
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(τὸ am) EY ERU) 
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= р -.: 


Fig. С-1 Input Impedance. 
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II. Radiation Resistance 
According to Chaney and Klamm [24], [25], the radiation resistance of a 


rhombic antenna in free space is given by 
К, = 120 {1.1544 + 2 Im (2.k1l xs а) + 2 Сї (2К1) - 2 С1(2К1 sin a) 

- Ci (2kl (1 + cos a) - Ci (21 (1 - cos a) ) 

2 ⁄ 

+ cos ( 2kl sin a ) сі ( 2k1 cos a (1 - cos a) )* Ci (21 sin a (l1 - sin a) ) 

T CI (2x1 cos a(l+cos a) ) = ci( ? kl sin a(l-sin a)) - 2ci ( 2 kl cows а ) 
Р m 
- 2ci ( 2kl sin а ) | 

+ ѕіп { 2к1 sin 8 Sil 2kl cos a (1-cos a) ) - Si ( 2k1l sin a (1-sin a), 
- Si (2x1 cos a (1+cos a) ) E SI ( 2x1 sin a (1+sin a) J + 2 Si C 2κ1 ee а ) 


T si ( at sin? a DED 


where 1 is the wire length - 2b - , a is the vertex angle oos WO 
Ci is the cosine integral 
x 


ci(x) = | = 


со 
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and Si is the sine integral 


x : 
Si(x) = Í ος dv 
0 


III. Input Impedance and Current Distribution 


Using equations (33) and dividing the line in sections as shown in 


Fig. D-l, we obtain for each section 


2. (να) ΠΤΙ Kave(I,n) B(l,n) X. Zin(b nd ОСТ) 


p Zin(I,n+1) C(I,n) + Kave(I,n) A(I,n) 


Each stub counts as only one section of uniform line. The number M 
of divisions within each radiation depends on the degree of approximation 
desired, but it is convenient to use for M an even number, to avoid the 


relatively faster rate of change of characteristics at the middle point. 
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APPENDIX D 


Computer Programming 
I. Program AZORES 

This program is written in code FORTRAN IV and it was prepared for 
the IBM 360 system that is operating at the Naval Postgraduate School. 

The program uses two subroutines. The main program reads the data, 
computes and writes the physical characteristics of the antenna, selects 
the current distribution for the antenna wires, controls the coordinates 
of the desired pattern, prints the results, and calls for the subroutines. 

The pattern computation is made by subroutine CORVO, that is given 
the antenna parameters, the elevation and azimuth angles, and the current 
distribution. In output, CORVO gives the two components of the electric 
field. 

The program computes each pattern twice, for two different models of 
current distribution. Model 1 (see Fig. 9) is fed in data cards; Model 2 
is generated by the main program. 

The use of a large number of comment cards, hopefully makes the 
program self-explanatory. As much as possible, the names given to 
algebraic and geometric symbols have the same phonetics as the ones used 
in the text. The only remarkable exception is the symbol TAU, that stands 
for p, the structure scaling factor. | 

The following pages contain the programming statements and the output 


print of a sample program. 
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II. Program SMIGUEL 


This program also is written in code FORTRAN IV for the IBM 360 
system. 

The program uses seven subroutines, as explained in SAMPLE PROGRAM 
SMIGUEL that follows. The subroutines supply to the main program the 
characteristics of the line and the non-uniformity corrections. The main 
program uses equations 33, first to find the input impedance, and then to 
find the current distribution. 

For each set of input data the impedance and current distribution are 
computed twice--one time without accounting for mutual impedance effects, 
and a second time accounting for them. 

Again, the program has a large number of comment cards that, hope- 


fully, make it self-explanatory. 
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